The Penrose-Fife phase-field model 
with coupled dynamic boundary conditions 

Alain Miranville* 
Laboratoire de Mathematiques et Applications, 
UMR CNRS 7348, Universite de Poitiers - SP2MI, 
Boulevard Marie et Pierre Curie, 
F-86962 Chasseneuil Futuroscope Cedex, France 
E-mail: Alain . Miranville@math . univ-poitiers . f r 

Elisabetta Rocca^ 
Dipartimento di Matematica, Universita di Milano, 
Via Saldini 50, 20133 Milano, Italy 
E-mail: elisabetta . rocca@unimi . it 

Giulio Schimperna* 
Dipartimento di Matematica, Universita di Pavia, 
Via Ferrata 1, 1-27100 Pavia, Italy 
E-mail: giusch04@unipv . it 

Antonio Segatti^ 
Dipartimento di Matematica, Universita di Pavia, 
Via Ferrata 1, 1-27100 Pavia, Italy 
E-mail: antonio.segatti@unipv.it 

January 24, 2013 



Abstract 

In this paper we derive, starting from the basic principles of thermodynamics, an extended version 
of the nonconserved Penrose-Fife phase transition model, in which dynamic boundary conditions 
are considered in order to take into account interactions with walls. Moreover, we study the well- 
posedness and the asymptotic behavior of the Cauchy problem for the PDE system associated to the 
model, allowing the phase configuration of the material to be described by a singular function. 
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1 Introduction 

In this paper we derive a model for phase transitions of Penrose-Fife type settled in a bounded domain 
O C K 3 . The peculiarity of our approach consists in the fact that we take into account the relations 
between fi and its exterior, including the effects of the interactions with the boundary into the free 
energy and entropy functionals ^ and S (cf. (f2"4"l) - (j2"5]) below) which drive the evolution of the system. 
A detailed derivation of the model is carried out in Section [2] The resulting PDEs system couples four 
nonlinear and singular evolution equations: two for the absolute temperature •& (one in the bulk f2 and 
the other on the boundary V) and two for the phase parameter X, which represents the local proportion 
of one of the two phases: 

f J. n, (i) 

|-A r (-I) + A' rW | + *(-i)^ro„r, (2) 
§ - AX - = in o, (3) 

^-A r X-V(X) + a,X = -^onr. (4) 

Here, A stands for the Laplacian with respect to the space variables in f2, Ar denotes the Laplace-Beltrami 
operator on T, and v is the unit outer normal vector to T. Moreover, A& and Ar are two quadratic functions 
of X related to the latent heat of the process, hb and hr are two heat sources, respectively in the bulk and 
on the boundary, and — s' Q b and — s' Q r are two nonlinear functions whose antidcrivatives — sqm and — so,r 
correspond to the configuration potentials of the phase variable. We admit the case when the domains of 
— so.b and — sq^ arc bounded, with the purpose of excluding the unphysical values of the variable X. In 
this case, we shall speak of singular potentials. A relevant example is the so-called logarithmic potential, 
given by 

- s(r) = (1 + r) log(l + r) + (1 - r) log(l - r) - -r 2 , S > 0. (5) 

With the choice ©, ±1 denote the pure states and the values X [—1, 1] are penalized by identically 
assigning the value — oo to s outside [—1,1]. In what follows, the potentials — s j, and — so.r will be 
split into the sum of (dominating) monotone parts /, and /r respectively, and of quadratic perturbations 
(cf. the last term in ([5])). The literature devoted to the mathematical features of phase transition models 
endowed with singular potentials is rather wide, also in specific relation with the Penrose-Fife model. 
Among the various contributions we quote [31] (for the Cahn-Hilliard equation), [21] (for the Caginalp 
phase-field model), [40] (for the Penrose-Fife system), and the references in these papers. 

A notable feature of system ([T])-6]) is the occurrence of dynamic boundary conditions. This type of 
conditions has been proposed in the literature in different contexts (for instance, in the framework of the 
Allen-Calm and Cahn-Hilliard models) with the aim of describing the interactions between the interior 
of a domain and the walls (cf., e.g. [13], [14], [15], [32], and references therein). In particular, the case of 
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singular potentials in the context of the Cahn-Hilliard evolution has been recently analyzed in [T7] , [T8] 
and [32], whereas the Caginalp phase-field system (cf. [6]) with dynamic boundary conditions on the phase 
parameter has been considered in several papers (cf., e.g. [7J, [5], [5], [IS]), while only recently it has been 
coupled with dynamic boundary conditions for both the phase parameter and the temperature. Such a 
problem has been considered in [50] and [TT], where the well-poscdness and the asymptotic behavior (in 
terms of attractors) of solutions have been studied, also in the case of singular potentials. 

As far as the Penrose-Fife model (cf. [34]) is concerned, a vast literature is devoted to the well- 
posedness (cf., e.g. [TO], [53], [35]) and to the long-time behavior of solutions both in term of attractors 
(cf., e.g. [25J, [35], [33]) and of convergence of single trajectories to stationary states (cf. [H]). Most of 
these contributions deal with Robin boundary conditions for the temperature, i.e., 

™ = f(*-*>r, (6) 

and no-flux conditions for the phase parameter: 

f =0on T. (7) 
av 

Relation ^ establishes that the heat flux through the boundary T is proportional to the difference 
between the internal and external temperature i? s > by a positive coefficient £, while ([7j prescribes 
that the boundary has no influence on the phase change process. Finally, we would like to quote the 
recent paper [26] where a phase separation model of Penrose-Fife type with Signorini boundary condition 
has been studied. 

The main novelties of our contribution stand in the fact that we can consider dynamic boundary 
conditions for a coupled phase-field system which may display a singular character both in i? and in X in 
view of the terms l/fi and s' or (X), s' b (X) in (jTJ)-(j4]) . Physically speaking, the main difference between 
the dynamic boundary conditions ([2]), (j4]) and the "standard" ones ©-([7]) stands in the fact that, in 
the case of dynamic b.c.'s, the walls of the the container have a significant effect on the phase transition 
process. This may include the case of two phase changing substances in contact, where one occupies a 
small layer surrounding the other one and this layer is approximated by a surface in the mathematical 
model (concentrated capacity phenomenon, see, e.g., [37], [38] for more details). 

A rigorous derivation of system flJ-Q starting from the basic laws of Thermodynamics is carried out 
in the first part of the paper. Namely, the equations are obtained by combining the free energy balance 
with the entropy dissipation inequality and imposing physically realistic constitutive expressions for the 
energy and entropy functionals. 

The second part of the paper is devoted to the mathematical analysis of the system in the framework 
of weak solutions. In comparison with similar models, the main mathematical difficulty consists here in 
the coupling between 

1. the singularity of the model, in particular due to the presence of 

(a) the term 1/d both in the bulk equations (fTJ) and ([3J and in the boundary equations ^ and (0J , 

(b) the (possibly) singular functions s' Q b in the phase equation ([3J in the bulk, and s' Q r in the 
phase equation ([4]) on the boundary, 

2. and the occurrence of dynamic boundary conditions. 
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Actually, as already noted in the case of the Caginalp model with dynamic boundary conditions (see 
|20j). if we are in presence of singular potentials, a particular care is needed. For instance, in order for 
equation © to make sense a.e. in the space-time domain, the monotone part / of — s' b needs to be 
controlled from above by the monotone part /r of -sj, r (cf. ((55)) below). In particular, this happens 
when the potentials — s' b and — s' o r have the same effective domain (e.g., the set [—1, 1], as in the case 
(|S])) and \s' ob \ explodes at most as fast as \s' or \ as X approaches the boundary of this domain (the values 
±1 in the specific case). 

It is worth noticing that, in the case of completely general singular potentials (i.e., without any 
compatibility condition assumed), a weak solution is still expected to exist; however, the equations ([3])- 

ruling the behavior of X need to be interpreted in a weaker sense either by means of duality arguments 
or of variational inequality techniques (cf. [2D], (32] for more details). This issue will be addressed in a 
forthcoming paper, where we also plan to weaken the regularity assumptions on the initial temperature 
$o- Actually, in addition to the natural conditions represented by the finiteness of the initial energy and 
entropy, we will assume here that i?o has the L 2 -rcgularity. This assumption, albeit meaningful, is not 
required in the physical derivation of the model and can be then considered to be somehow artificial. In 
other similar contexts (cf., e.g. [30]), the L 2 -rcgularity of the initial temperature has been shown to be 
not necessary and, for this reason, we will try to remove it also for this model. 

After proving well-posedncss of system ([T])-((4]), we shall analyze further properties of solutions. Pro- 
ceeding along the lines of [3D] where the case of standard b.c.'s is treated, we shall prove, by using a Moscr 
iteration argument, that the temperature is uniformly bounded from below for strictly positive times. 
Moreover, if -do is slightly more summable (cf. below), we also have a uniform upper bound. This 

kind of behavior occurs commonly in the framework of parabolic equations with very-fast diffusion terms 
(see [42] and [3] for the Cauchy problem in the whole space and the recent contributions [40] and [41] for 
the bounded domain case with Neumann and dynamic boundary conditions, respectively). In particular, 
in dimension three the exponent p = 3 appears to be critical in the sense that solutions starting with 
initial conditions in L p with p > 3 become L°° for strictly positive times. On the other hand, for p < 3 
the situation is drastically different as the self-similar solution in R 3 ((•)+ denoting the positive part) 

shows. The smoothing effect for p = 3 is currently an open problem (see [42] for further details). Moreover, 
it is worth noting that, at least when no external source is present, the regularization estimates are also 
uniform with respect to time and give rise to additional regularity properties for both components of the 
solution. 

Taking advantage of the regularization estimates, in the last part of the paper we finally investigate 
the long-time behavior of trajectories. Namely, we are able to prove that, at least in case (|111[) holds, 
any weak solution admits a non-empty w-Umit set which only contains stationary solutions of the system. 

A more precise characterization of the long-time behavior relies on the structure of the set of steady 
state solutions, which requires some further explanation. Actually, integrating ([T]) in space, using the 
boundary condition ([2]), and assuming zero external source, one readily sees that the value 

H= I ($ + \ b (X))dx+ [ (0 + Ar(X)) dS, (9) 
Jo. Jt 

representing the "total mass" of the internal energy, is conserved in the evolution of the system. Conse- 
quently, any limit point of a given weak solution has also to respect the constraint ([9j), with /i depending 
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only on the initial datum. However, when /i is small, we are not able to exclude that the set of stationary 
states satisfying ([9]) might contain temperatures $ being arbitrarily close to (note that the stationary 
formulation of ([T])-([2]) simply prescribes that $ is a constant function). Consequently, the right-hand sides 
of the steady state system associated to ©-(EH might be arbitrarily large, which considerably weakens 
the regularity properties of the set of its solutions. 

On the contrary, it is easy to show that this situation cannot occur when either \x is large enough or 
A' satisfies a suitable sign condition (cf. (|214l) below). If cither property holds, the set of stationary states 
is bounded in a very strong norm, independently of the magnitude of the initial data (but depending on 
the value of fj,). This fact, together with the precompactness of solution trajectories in the natural energy 
space and with the existence of a coercive Liapounov functional (namely, the energy ty), implies that the 
system (HJ)- admits a smooth global attractor, which is the last result we prove. 

We conclude by giving the plan of the paper: in the next Section[2]we detail the physical derivation of 
the model from the basic thermodynamical principles. In Section [3] we introduce our precise concept of 
solutions and state our main mathematical results related to well-poscdness and regularization properties 
of weak solutions. The proofs are given in the subsequent Section 2J Finally, the long-time behavior of 
solutions is separately analyzed in the last Section [3] 

2 Derivation of the model 

The Penrose-Fife model is derived by considering the free energy density w and the entropy density s, 
assuming that these quantities depend both on the order parameter X and the absolute temperature d 
(as in the original paper by Penrose and Fife (34]), and imposing that the basic laws of Thermodynamics 
are satisfied. Assuming that we always have sufficient regularity of the involved variables, we can then 
write 



w = e — i9s (Gibbs' relation), 



(10) 



where 



s 



dw 



(11) 



and the internal energy density e is defined by 



e 



(12) 



Moreover, we assume that the total free energy has the expression 




(13) 



whereas the entropy functional is given by 




(14) 



where f2 is the domain occupied by the system and k > denotes an intcrfacial energy coefficient. 
The evolution equations for X and d are then obtained by stating the relations 



% = KAf e+ U,K>0, (16) 

where 5 denotes a variational derivative and ft, is a source term. 
Taking finally 

e = c i9 + A(X), (17) 

where Co > stands for the specific heat of the system (assumed constant) and A is the latent heat 
density, typically (cf. [34] p. 53]) given by 

A(r) = ~ar 2 + br + c, a > 0, (18) 

we find (see, e.g., [30]; see also below) 

s(e,X) = colntf + s (X) +a. (19) 

Here, —sq denotes a configuration potential, typically having a double-well character (one can also con- 
sider a logarithmic double-well potential of the form (O), and such that Sq < S, S > 0, and c\ is a 
constant. These choices give rise to the Penrose-Fife system: 

f = K .(,AX + 4W -^), (20) 

These equations are usually endowed with the boundary conditions (cf. the Introduction) 

dX 

— = on T, (22) 

|* = _ f (0 - a ) on T, £, <? s > 0, (23) 

where T = <9f2 and ^ is the unit outer normal vector to Y. 

Now, in order to take into account the interactions with the exterior of f2 (e.g., the walls), it is natural, 
following [27] (see also [IT] and [H]), to add a boundary contribution to the total free energy and to take, 
in place of (fT3j) . 

* = *(X,VX,V r X,?9) = ^(^rl VX ! 2 +ui b (4,X)j dx (24) 

(^|VrX| 2 + wr(^,X)^ dS, k 6 , *r>0, 

where Vr is the surface gradient and and wr are the bulk and surface free energy densities, respectively. 
Similarly, it is reasonable, in view of ([24]) . to introduce the generalized entropy functional 



S = S(X, VX, VrX, e) = J (-y |VX| 2 + s b (e, X)) da- + J (-^ |V r X| 2 + s r (e, X)) d£, (25) 
where and sr are the bulk and surface entropy densities, respectively. 



G 



As before, we assume that Gibbs' relation holds, namely, 

Wb = e — #s& in 0, 



and that 



sr 



$sr on r, 



dwr 



in i7, 

on r, 

in fl, 

on r. 



We now note that, in view of (125 



SS . ds r dX 
- = KrAr X + — - Kb - 



on r, 



where Ap is the Laplace-Beltrami operator. 

Then, assuming that, as in the classical model, 



dX 
~dt 



K* 



5S_ 

sx' 



K* > 0, 



we obtain the equations 



OX 



dt 



ds b 



— = K* [ n b AX + — - 



ox 



dx 
~dt 



K* n r A r X 



9sr 



Kb- 



in fl. 

dX 



on r. 



(26) 
(27) 

(28) 
(29) 
(30) 

(31) 

(32) 
(33) 

(34) 

(35) 
(36) 



dX u dv 

Next, in order to describe the evolution of the temperature, we generalize (fTB"]) as follows. We intro 
'win 



duce, for U 



(regular enough at this stage), the linear operator A defined by 
AU = 



-Au|n 
-A rW | r + ^| r 



and write that 







en 



dt \e r 



= —KA 



SS_\ 

Se \ n 
SS I 



H, 



(37) 



(38) 



where H 



is a forcing term, with h b and hr standing for the bulk and boundary heat sources, 



respectively. Introducing the linear operator A is natural when considering dynamic boundary conditions 
(see [IT] and HU). Noting that 

— = -T- m S2, (39) 
oe oe 
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Se de 

we deduce the equations 

— = A A— + h b m 11, 
at oe 

de . dsr , T ^ d ds b ^ 
— = KA r —^- + h r -K — — - onT. 
at oe ov oe 

We then assume that (see (|17p) 

e = c , 6^ + A 6 (X) in Q,, c ,6 > 0, 

e = c<),r$ + Ar(X) on T, c ,r > 0. 
Noting that it follows from ([3D]) and ([3T]) that 

e = -^onr, 
we have, freezing X and integrating between i?o and #, &q, $ > 0, 

w 6 (tf,X) w b {$ ,X) [* r Co, b \ b (X) . 

q 5 = _ / ( 1 2~ ) dT m S2 ' 

t«r(<?,X) w r (^ ,X) f%.r,M4, r 
S q = -/ + — dronT, 

which we rewrite in the form 

w b (0,X) = -co, 6 tflnt? + c , 6 tflntfo + A 6 (X) - fo , b {X) in Q, 

ior(i9, X) = -co,ri? lm? + co,ri? In i? + A r (X) - i?s ,r (X) on T, 

where we have defined s ,6 := (wh^OjX) — Af,(X))/#o and s ,r := (wr($0jX) — Ar(X))/#o- 
We finally deduce from (|28 > pg ]) and (|3g ]) -([50 ]l that 

Sfe(e,X) = c ,&ln# + s ,&(X) + ci, fc , ci,& = -c ,hlntf + c ,fc in Ct, 

s r (e,X) = c , r hi^ + s ,r(X) + Ci, r , Ci, r = -Co.rln^o + co,r on I\ 
We now note that it follows from ([IS |) -(|33 p and (f5T ]) -(f52 |) that 



• e- A 6 (X) ' 

CQ, b 

a r (e, X) = c .r In ( 6 ~ Ar ^ ) + s ,r(X) + c 1T on T, 



s&(e, X) = co.h In f ) + so,&(X) + ci,6 in O, 

\ co, b ' 

■e-x T {xy 

cot 



which yields 

+ s' 0tb (X) in fi, (55) 

+ s' 0iV (X) on T. (56) 

We finally deduce from ([3U)) - (|4"2")) and ([55j) - ([5"6"| the following Penrose-Fife system with dynamic bound- 
ary conditions: 





1 


ds b 


Km 


5e 


~ ' 


dX ~ 


■d 


9sr 


1 


dsr 


X' r (X) 


<9e 


~ 0' 


~d~X = 





— = A'* U b AX + s'^ b (X) - in 0, (57) 

9X - ^* (^rArX + 4 r (X) - - k 6 |^ on r, (58) 



at " y*— 1 -- 1 -o,rw tf 

c , b ^ = —KA— - A' 6 (X)— + ft 6 in 0, (59) 

co, rii - = -KA r - - X' r (X)— + h r + K-- on V. 60 
at v at av v 

Taking, for simplicity, all constants equal to one, (|57[) - (|6"0"|) reduces to 

in n, (61) 



AX + s' b {X)-^^h 



at u > ' 

dX i , . XUX) dX , . 

-=A r X + s' oA X)-^>--ouT, (62) 

Brf i Fix 

- = -A--X b (X)— + h b mQ, (63) 



av . i w , .ax , a i 

which can also be rewritten in the following compact form 



A r-- A rM— + ^r + — - on T, (64) 



sv^r^UW + ^wJ-^J' (65) 

The remainder of the paper is devoted to the mathematical analysis of system (|6"5|) - (f6"6")) in the framework 
of weak solutions. 
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3 Main assumptions and preliminary results 



We introduce here our main assumptions, together with several mathematical tools which are needed in 
order to provide a precise analytical statement of our results. 

We let f2 be a sufficiently smooth, bounded, and connected domain in R 3 with boundary T. We set 
ft := il U r, set H := L 2 (il), and denote by (•, •) the scalar product both in H and in H 3 and by || ■ || 
the related norm. Next, we set V := i? 1 (fi) and denote by V the (topological) dual of V. The duality 
between V and V will be indicated by (•, •}. Identifying H with H 1 through the scalar product of H, it 
is then well known that V C H C V with continuous and dense inclusions. In other words, (V, H, V) 
constitutes a Hilbert triplet (see, e.g., [29]). Such a triplet is usually used for stating weak formulations 
of elliptic or parabolic problems defined on f2. 

However, since system (fG5 ]) - ([66"l) also includes equations defined on T, we need to introduce some 
further spaces taking also boundary contributions into account. Thus, we set Hr := L 2 (T), Vr := 77 1 (r), 
and denote by (•, -)r the scalar product in Hr, by || • ||r the corresponding norm, and by (•, -)r the duality 
between Vf and Vr. In general, the symbol || ■ \\x indicates the norm in the generic (real) Banach space 
X and (-,-)x stands for the duality between X' and X. We also denote by Vr the tangential gradient 
on r and by Ar the Laplacc-Bcltrami operator. Thus, we can define the spaces 

H:=HxH v and V := {z G V : z\r £ Vr}- (67) 

Here and in the following, z\r, or also zr, will denote the trace of z in the sense of a suitable trace 
theorem. Next, we introduce the 7^-scalar product in the following natural way: 



((&, k), (s, a)) := (k, s) + (k, a) 



r ■ 



Then, we set, respectively on V and on V , 



j (Vz ■Vffl) da; + / (z\ T w\r + V r ^|r • V r w|r) dS, (69) 
Jn Jr 

((z,w)) v := { {Vz-Vw)dx+ [ z\ r w\ r dS, (70) 



[z, w)) v := 



together with the related norms || • ||y, || • It is not difficult to prove (see, e.g., [33l Lemma 2.1]) that 
the space V is dense in %. Concerning the scalar product in (|70|) . we notice that the related norm || • ||y 
is obviously equivalent to the usual one. We also set 

W := {z e V : ze H 2 (n), z\r € H 2 (T)} (71) 

and endow this space with the graph norm, so that W C V, continuously and compactly. 
The above defined functional spaces allow to introduce some elliptic operators. We set: 



A:V^V, (Azi,z 2 ) := / Vzi ■ Vz 2 dx, (72) 

A r : V r -+ V£, (A r Cu C 2 )r := J V r Ci ■ V r C 2 dS, (73) 
A-.V-+V, (Az 1 ,z 2 )v ■= (A Zl ,z 2 ) + (A r zi,T,Z2,r)r- (74) 



We shall also use the operator A defined in (|37|) . Note that A can be interpreted as an operator defined 
on VV and taking values in "H, thanks to the trace theorem for normal derivatives. 
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In what follows, we will use the following convention: as far as equations on O are concerned, the 
elements of V will be interpreted as functions defined on il with the proper regularity. When, instead, 
as in most cases in the paper, a system defined on f2 x T is considered, then the elements of V will be 
considered as pairs of functions (z, z|r)- In other words, V will be identified with a (closed) subspace of 
the product space iJ 1 (17) x -ff 1 (r). Analogously, V will be identified with a subspace of H 1 (f2) x i/ 1 / 2 (r), 
in view of the trace theorem. If we have, instead, h £ H, h will be often thought as a pair of functions 
belonging, respectively, to H and to Hr, and both denoted by the same letter h. Of course, if we do not 
have additional regularity, the second component of h needs not be the trace of the first one. Identifying 
H with H' through the scalar product we obtain the chain of continuous and dense (thanks to 

the density of V into H, to the density of H 2 (£l) into V, and to the continuous inclusion H 2 (£l) C V) 
embeddings 

VcVcHcV'cV. (75) 

In particular, the relation 

((k,K),(z,z\r)) H = / kzdx+ / nz\ r dS = ((k, «), z)v (76) 
Jn Jr 

holds for any z £ V and (k, k) £ H. Of course, an analogous relation could be stated for z £ V. 
Next, for any function, or functional z, defined on Q, we set 



(77) 



where the integral is substituted with the duality (z, 1) in case, e.g., z £ V . We also define the measure 
dm given by 

/ vdm:= / vdx+ / Vr dS, (78) 
JTi Jn Jr 

where v represents a generic function in L 1 ^) x X 1 (r). With some abuse of notation, we will also write 

m{v) := W+W\ ii udTO ' (79) 

i.e., the "mean value" of v w.r.t. the measure dm. Here |T| represents the surface measure of T. 

With these functional spaces at our disposal, we can now state our hypotheses on the nonlinear terms. 

For convenience, we split — s' b (respectively, —s' r ) into a sum of a (dominating) monotone part / 
(respectively, /r) and a linear perturbation. More precisely, we assume that 

- s' , b = f(r) - 5r Vredom/, -^ >r (r) = /r(r) - Sr Vr £ dom/r (80) 

and for some 5 > 0, with 

/ £ C°(dom/,R), / r G C° (dom/r, M) monotone, /(0) = / r (0) = 0, (81) 

where dom/ and dom/r (i.e., the domains of / and /r) are open intervals of R containing 0. We will 
say that, for instance, — so,fc is a "singular" potential if its domain does not coincide with the whole real 
line (and we will say that it is a "regular" potential otherwise). In both cases, we will assume that 

lim (/(r) — Sr) signr = lim (fr(f) — Sr) signr = +co. (82) 

r— >c? dom / r~^d dom /r 
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The key assumption that will allow us to obtain a pointwise estimate of the nonlinear terms /(X) and 
/r(X) is the following compatibility condition: we ask that there exist two constants c s > and C s > 
such that 

dom(/ r ) C dom(/), f(r)f r (r) > c s \f(r)\ 2 - C s Vr G dom(/ r ). (83) 

In other words, the boundary nonlinear term fr(f), up to the sign, has to be larger than the bulk nonlinear 
term f(r), at least for r far from 0. We also introduce, whenever they make sense, the antiderivatives 

F(r) := f f(s) ds, F r (r) := f f T (s) ds. (84) 
Jo Jo 

Notice that, in case (for instance) dom/ is bounded, but / is globally summablc on dom/, F can 
(and will) be extended by continuity to dom/. This is the case, e.g., of the logarithmic potential ([5]). 
Moreover, F (and, analogously, Fp) will be thought to be further extended to the whole real line by 
assuming the value +oo outside its effective domain. Then, identifying / and /r with maximal monotone 
graphs in R x K, we have / = dF and fr = dFp, d representing the subdifferential of convex analysis 
(here in K.) . We refer to the monographs [TJ [21 H] for an extensive presentation of the theory of maximal 
monotone operators and of subdiffcrentials. 

To obtain an estimate of the full I^-norm of u, we will also need a proper form of Poincare's inequality 
(see, e.g., [40j Lemma 3.2]): 

Lemma 3.1. Assume that il is a bounded open subset ofM. d . Suppose v G W 1 ' 1 ^) and v > a.e. in 
il. Then, setting K := J„(log«) + dx, the following estimate holds: 

\\v\\ LHa) <\n\e c * K + ^\\Vv\\ L i (a) , (85) 
the constants C\ and C2 depending only on fi. 

Besides assumptions (|80j) - (f83j) , wc shall analyze system (|6H - (|64j) under the following hypotheses: 

\ b E C 2 (dom/) with Aft G L°°(dom/), A r G C 2 (dom/ r ) with Ap G L 00 (dom/ r ), (86) 

(i?o,r?o) G H, ^o,Vo > a.e., (log^o, log?7o) G L^Sl) x i 1 (r), (87) 

X G V with s 0:6 (X ) G L 1 ^) and s 0: r(^or) 6 ^(r), (88) 

H G L 2 (0, T; H), m(H) = a.e. in (0, T). (89) 

Let us define the energy functional as 

S[(d,Ti),(X,Xr)] ~J(#-log# + \ b (X) + ^--s 0tb (xj) dx 

+ J(v~ log' 7 ? + A r(X r ) + |Vr ^ r|2 - flo,r(X r )) dS, (90) 

whenever it makes sense. We will often just write f [t?, X] for brevity. Note that we admit £ to take the 
value +00. However, we can readily check that the above assumptions ([57| - ([55)) make the initial energy 
finite. Namely, we have 

E := £[(0 O , Vo ), (X , X or )] < 00. (91) 
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Together with (j8"o) - (|59")) . we will also need the energy functional £ to be coercive with respect to (X, Xr). 
This is obtained asking that there exist cy > and c £ K such that 

£[&,X] > c v \\X\\l -c for all pairs (tf,X). (92) 

Such an assumption is satisfied for instance whenever 

Afc(r) — so,b(r) > c\r 2 — C2 for all r € dom/ and Ar(V) — So,r(?") > ciT 2 — C2 for all r <E dom/r, (93) 

and for some c\ > and C2 € K. Of course, (|92[) holds trivially if — so.t, or — So,ri or both, are singular 
potentials, or also in the case of standard double-well potentials. 

It is worth noting that the first of assumptions (|87[) on ($0, Vo) is somehow artificial, in the sense that 
it is stronger than what would be required in order for the initial energy to be finite. Actually, one could 
relax (f57|) by taking just (tfoi^o) G i 1 (0) x i 1 (r) or even (#o,?7o) G V, paying the price of having a 
weaker (and more delicate to deal with) notion of solution. The Neumann-Neumann Penrose-Fifc model 
with L 1 or V' initial temperature has been recently studied in |40j . In a forthcoming paper, we intend 
to analyze also the present model in a similar regularity setting. 



3.1 Weak solutions 

We use the functional framework introduced above to specify a rigorous concept of weak solution to the 
initial value problem for (|65[) - (|66l) . named Problem (P) in what follows. 

Definition 3.2. We say that a quadruplet (■&, n, u, X) is a weak solution to Problem (P) if there hold 

(0,-n) e H\Q,T; V) nL°°(0,T;H), (94) 

(logt?, logr?) G L°°(0, T; L\Q) x L 1 ^)), tf,n>0 almost everywhere, (95) 

u e L 2 {0,T;V), (96) 

X e i? 1 (0,T;-H)nL oo (0,T;V)n J L 2 (0,r;W), (97) 

/(X) 6 L 2 (0, T; H), /r(X r ) G L 2 (0, T; H r ), (98) 

together with, a.e. in (0,T), the equations 

i i 

u = — — a.e. in D,, n = a.e. on T, (100) 

v ur 

d f x \ , a{ x \ ( s ibW \_f A' b (X)u 



and the initial conditions 

(0,»?)|t=o = (0o,»K>), (X,X r )|t=o = (Ao,Xo r ) a.e. in and on T. (102) 
Sometimes, for brevity, we shall indicate a solution just as a pair (■&, X) rather than as a quadruplet 

(tf,T?,U,X). 
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Remark 3.3. It is worth giving some explanation on the boundary behavior of i9. Since u = — l/i9 £ 
L 2 (0,T; V) by (f96|) . it turns out that u has a trace ur on T, which belongs to Vr for almost every value 
of the time variable thanks to the definition of V. On the other hand, we cannot simply write rj = ??r 
since the trace of i? needs not exist. We have, instead, to intend rj as (minus) the reciprocal of the trace 
of u, as specified by the second (jlOOp . When we consider smoother solutions (for instance, in the a priori 
estimates, or in the case when #q is more summable, cf. below), this problem does not occur since 

the higher regularity of i? permits to give sense to its trace. Regarding the phase variable, the situation 
is simpler; indeed, by ([9"7| . X is always smooth enough to have a trace Xr- 

3.2 Existence and uniqueness results 

In this part we state our main results regarding well posedness of Problem (P) and regularization proper- 
ties of weak solutions. In what follows we will denote by c a positive constant, which may vary from line 
to line (or even in the same formula), depending only on the data of the problem. Specific dependences 
will be indicated when needed. Moreover, we will denote by Q a nonnegative-valued, continuous and 
monotone increasing function of its arguments. Our main result can be stated as follows: 

Theorem 3.4. Under assumptions ([50)1 - ([53]) . ([5B"]) - ([59"[) , and ([9"2"| . there exists a unique weak solution to 
Problem (P). Moreover, for all t > there holds the energy identity 

£[(0(t),»7(t)), (X(t), X r (t))] + f (l|V«(a)|& + \\X t {s)f n ) As = E + f (H(s), u(s)) n da. (103) 

Jo Jo 

If, in addition to ([89|h the forcing function satisfies 

H = {hb, h r ) G i? X (0, T; V') n L 2 (0, T, L 3+£ (ft) x L 3+e (r)), (104) 
then we also have, for any r G (0, T), the following regularization properties: 

\\u\\ L °°(t,T;V) + ||«||x,«((T,T)xn) < 0(Eo, r' 1 ), (105) 
l|X|Uo= (TiT .^ (n)) + \\f{X)\\ L oo (TtT . H) < Q(Eq, r" 1 ), (106) 

hr\\ L °°(T,T;v T ) + ll«r||i-((r,r)xr) < Q(Eo, r- 1 ), (107) 
||Xr|U=o (TiT;H2(r)) + ||/r(X r )||Loo (TiT;ffr) < Q(E , r" 1 ). (108) 

The proof of the existence part is based on an approximation-a priori bounds-passage to the limit pro- 
cedure. Additional estimates lead to (|105[) - ([108[) . In particular the L°°-bounds in (|105p and in (|107|) 
are obtained by adapting a Moser iteration scheme with regularization devised in [ID]. Note that the 
existence part generalizes to the dynamic boundary conditions case the result of Ito and Kenmochi [24] . 

In the case of singular potentials, we can also prove that, at least for strictly positive times, X is 
uniformly separated from the singularities of the potentials. In order to avoid unnecessary technical 
complications, we shall state this property under the additional assumption that 

dom/ = dom/ r = (-l,l), |/ r (r)| > «,|/(r)| - C s Vr G (-1, 1) (109) 

and for some k s £ (0, 1], C s > 0. Namely, we assume the potentials to be normalized so that the pure 
states are represented by the values ±1, both on the bulk and on the boundary. Moreover, in view of 
(|83p. we require |/r| to be larger than (a positive constant times) |/|, at least in proximity of ±1. 
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Corollary 3.5. Let the assumptions of Theorem 13.41 hold, together with (|109p . Then, for any r > 0, 
there exists lo £ (0, 1), depending on t but independent ofT, such that 

\X(t, x)\<l — w almost everywhere in (r, T) x Q and in (r, T) x T. (110) 

Next, we discuss the asymptotic regularization of the temperature field. Actually, we have already noted 
(cf. (|105|) . (HU7|)) that condition (|104p on the forcing function H is sufficient in order for $ to become 
uniformly separated from zero for strictly positive times. The following result (which generalizes |40[ 
Thm. 2.7], where Neumann conditions are considered) states that d is bounded from above, at least 
for strictly positive times, provided that (|104[) holds and the initial temperature enjoys some additional 
summability property. 

Proposition 3.6. Let the assumptions of Theorem l3.4l hold. Let also assume that 

(0 O) Vo) G L 3+e (n) x L 3+£ (T), for some e > 0, (111) 
cither (fTM|) holds, or X £ L°°{Ti,dm). (112) 

Then, any weak solution to Problem (P) satisfies, for any t £ (0,T), 

||tf|U°=((T,r)xf2) + \\v\\L°°((r,T)xr) < Q(Eo,r _1 , l^ollz^r), IMU*+«(r)i l|Xo||£=o(n )dm) ). ( 113 ) 

We point out that, in the case of singular potentials, the uniform boundedness of Xo required by (|113j) is 
a direct consequence of 

As a byproduct of Prop. 13.61 we can also prove additional regularity of the time derivative of 
Corollary 3.7. Under the assumptions of Proposition 13.61 we have 

$£H 1 (T,T;U)foranyT>Q. (114) 
As a consequence, (|99|) can be decoupled and interpreted in the strong form (I66|) as a relation in TL. 

4 Proofs of the main results 

4.1 Proof of Theorem 13. 4t a priori estimates 

As a first step, we detail the main estimates constituting the core of the existence proof. In order 
to simplify the exposition, we limit ourselves to perform formal a priori bounds on the solutions of 
Problem (P). In the next section we will see that these bounds imply weak sequential stability. It is clear 
that, in a formal proof, the estimates should be performed in the framework of a proper approximation 
scheme (e.g., a Faedo-Galerkin approximation or a time discretization), possibly combined with some 
regularization of the data. However, this kind of procedure has been already described in full detail in 
several papers related to similar models (see, e.g., |17j and |20j ) and, actually, the arguments given in 
these papers could be easily adapted to our case. 

That said, we start by presenting the estimates. In all what follows, we shall assume to have sufficient 
regularity to justify all the computations. In particular, we ask ■& to be smooth enough to have a trace, 
so that rj = i?r (and correspondingly 1/rj = l/i?r)- 
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Energy estimate. Test (|99|) by ( J I/77J an< ^ 051]) by ( ^ ' )■ Noting that two terms cancel out 



J ~ J/Jj and I 
and using that H has zero mean value, we obtain 

+ J |V(-l/i?)| 2 d.T + J^\V r (-l/r,)\ 2 dS + JjX t \ 2 dx + J^\X Fit \ 2 dS = -(H, (115) 

Using that m(H) = (cf. ([75])), we can write 

(H, = (H, 1/0 - mn(l/tf)) w = j n hb {~ mn Q + l) dx + J r ' lr {~ m °G) + J) d5 - 

Now, the integral over f2 is easily estimated, using the Poincarc-Wirtinger inequality, as 

^(i-^nQ) cL^r < c e ||^ b || 2 +e||V(l/^)|| 2 Ve>0. 
On the other hand, we treat the integral over T in this way: 

J fe r(^-m fi Q) dS< ||fcr||r||l/i7-mn(l/0)||r 

< c||/i r ||r||l/^ - mn(l/tf)|| v < c||fc r ||r(||l/tf - mn(l/0)|| + ||V(1/0)||) 

< c||Mr||V(l/tf)|| < c e ||Mr + e||V(l/0)|| 2 Ve > 0, (116) 

where in the second and in the fourth inequalities we have used, respectively, the trace theorem and the 
Poincarc-Wirtinger inequality. Hence, taking e small enough and integrating (|115p . we obtain 

£{t) + ( ||V(1/,9)I[2 + H^CV^IIr + nM\u) ds < E + c||H||2 2(0>T . w) . (117) 

To control the full U-norm of l/fl, we have to provide a bound of its mean value (actually, only the 
gradient is estimated in (|117[1 ). To this aim, we use Lemma 13-11 with v = — l/i?, obtaining 

\\l/n LHn) < \Q\e c ^^r + ^||v(l/tf)|| il(n)> (118) 

and the first term on the right-hand side is uniformly bounded thanks to (|117[) and to the expression (|90j) 
of the energy functional. Moreover, by the trace theorem, we get 

IIV»/llvr<c(||l/i?||v + ||Vr(l/»7)||r). (119) 
Hence, collecting the above computations, (|115[) gives the a priori bound 

S[0,X](t)+ [ T (\\(l/#)(s)fv+\\Xt(s)\\K) d S <g(E ,T,||H||2 ) Vt<T. (120) 
Jo 

Estimate of the nonlinear terms. Estimate f|120[) gives that u = —l/fl E L 2 (0,T;V) and X £ 

L°°(0,T;V), implying, thanks to (|85|) . that the right-hand side of the phase equation (|10ip belongs 



16 



to L 2 (0,T;H). Now, let us test (fTUTJ) (both on the bulk and on the boundary) by f(X). Using the 
monotonicity of / and assumption (|86j) we then get 

j f J_F(X) dm+||/(X)|| 2 + Jj(x r )(f r (X r ) - 5X r + dS 

< jf /(X)(«5X - dx < i||/(X)f + c||X|| 2 + c(l + ||X||£< (n) )Hli« (n) . (121) 

The key point is represented by the control of the last term on the left-hand side, and here the compati- 
bility condition (|83[) comes into play. Indeed, using (|83[) . we readily arrive at 

J t f(Xr) (/r(Xr) - SX r + > |||/(X r )|£ - c c||X r || 2 - c(l + ||X r || 2 L4(r) ) Kll^(r)- (122) 

Hence, integrating (|121[) in time and using (|122[) and f|l 20[) . we infer 

\\f(X)h*(o,T;H) + ||/(Xr)||L»(o,riJrr) < Q(Eo,T). (123) 

Here and below, we allow Q to depend additionally on H. Moreover, we recall that ll-FX^o)!^^ dm ) < 
Q(E ) thanks to (j88|) and (f83|) . 

Next, (|120|) . fjl 23(1 . and a comparison of terms in (foTj) (i.e., the bulk component of (|99|) in the strong 
formulation - recall that we assume the solutions to be smooth at this level) , give 

||AX|| i2(0iT . H ) <Q(Eo,T). (124) 

Consequently, using standard trace and elliptic regularity theorems (cf., e.g. [5] Theorem 2.7.7 and 
Theorem 3.1.5]), it is not difficult to arrive at 

\\d„X\\ L 2 (0!T . Hr) <Q(E ,T). (125) 

This allows to test the boundary equation (p2"|) by /r(Xr). Proceeding as above (but without using (|6T|) ) 
and controlling the term d v X directly by means of (|125p . we finally arrive at 

l|ArXr|| L 2 (0!T . Hr) + ||/r(X r )|U»(o,T;ir r ) ^ Q( E o,T)- (126) 

Now, let us come to the temperature equation. By (|120j) and assumption (|86|) . we have 

\\dth(X)\\L2(o,T;L3/i(n)) < ll^( x )IU~(o,T;L6(o))||X t || i 2( ^ T;fl - ) < Q(E ,T), (127) 

and a similar relation on T. Actually, when — so,6 an d — So,r arc singular potentials, we automatically 
have a uniform L °° -bound on X. Hence, we obtain more precisely 

\\dt\b(X)\\mo,T;H) < ||A , 6 (X)|| i oc (Q) ||X t || i 2 ( o,T;H) < Q(E ,T). (128) 

By a comparison of terms in the (coupled) weak formulation (|99]l. we then get in both cases 

IMU»(o,t;VO <<9(Eo,T). (129) 
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Filially, to get the ^-regularity of tf, we test by •&. Then, using the first of ([57 ]) , ([55 ]) and ([125]) . 
noting that 



< c(l + ||X|| L< » (n) )||X t ||||^|| < c||X t || 2 + c(l + ||X|& a(n) )||tf|| 2 , (130) 



thanks also to (|86p , observing that a similar relation holds on T, and applying Gronwall's lemma, it is 
not difficult to arrive at 

\\n L oo {0 , T - m <Q(Eo,T). (131) 
Relations (jT^ and (|T3"Tj) give the desired (jMj) . 

4.2 Proof of Theorem 13. 4t weak sequential stability 

In this part, we prove weak sequential stability for Problem (P). Namely, we assume to have a sequence 
(i? n ,X„) of sufficiently smooth solutions having uniformly bounded initial energies Eo, ra and show that, 
at least up to the extraction of subsequences, this family converges in a suitable sense to a weak solution 
to the problem. Of course, in principle (i? n ,X„) could also be taken as the solution to some regularized 
Problem (P ra ); in this respect, the procedure below can be seen as a simplified version of the arguments 
needed to remove the regularization. However, as noticed at the beginning of Subsec. 14.11 we decided, 
for the sake of simplicity, to skip the details of the approximation procedure. 

That said, set u„ := — l/# n . Consequently, its trace M nj r is equal to —l/r) n . Thanks to estimates 
(|120|) . (|124j) . (|126|) and (|129|) - ()131[) . an d to the assumed uniform boundedness of the initial energies, we 
then obtain 

sup {||w n || L 2 (0 iT;V ) + ||9 t X n j| i 2( 0:T; - H) + ||X„|| i 2( T;W) + ||so 6 (X n )|| i 2( 0: r ; ff) (132) 



+ ll s 0,r(Xri,r)||L 2 (0,T;ff r ) 



where c is a positive constant independent of n. Then, using standard weak compactness arguments and 
the Aubin-Lions lemma, there exist u, X with traces, respectively, ur, Xr, and a non- relabelled subsequence 
of n such that 



U u weakly in 1/(0, T; V), (133) 



X n X weakly in H l {Q,T-H) f)L 2 (0,T;W) and strongly in L 2 (0, T; V), (134) 

»o,&(*») i ^ ^7 weakly in L 2 (0, T; if), (135) 



o,r(X„) ^-^> weakly in L 2 {0,T;H r ), (136) 



Ai(X„)«„ ^> V 6 weakly in L 2 (0,T;H), (137) 
Ar(X„ >r )M ra ,r % weakly in L 2 (0, T\ H v ). (138) 

Then, thanks to the properties of A and Ar, it is not difficult to obtain that A'(X n ) — > A'(X) strongly 
in L 2 (0,T; H) (similarly for Ap). Moreover, by the usual monotonicity argument [2j Prop. 1.1, p. 42] 
(recall that s' b and s' r verify ([80]) ). we can identify, almost everywhere in (0, T) x fl (or in (0, T) x T), 
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s' ob = s' ob (X), s' Q r = s or (Xr), X' b = X' b (X)u and A r = A r (Xr)wr- Now, we can take the limit in ([9"9")l . 
Indeed, by (|131|) we have 

(0n,0n,r) weak star in L°°(Q,T;H). (139) 

To conclude we have to show that — l/i? = it a.c. in (0,T) x SI and that — 1/77 = u-p a.c. on (0,T) x I\ 
Actually, these identifications are consequences of the monotonicity of the map v *-> —1/v (and of its 
realization in L 2 (0,T;H)). Indeed, by (p3]l , 

m„ ^^^^ u weakly in L 2 (0,T; V) and consequently in L 2 (0, T; L 6 (S1) x L 6 (r)), (140) 



Moreover, noting that, thanks to (|129[) and (|131[) . 

tf„ strongly in L 2 (0, T; V), (141) 

we get the following limsup-incquality: 

limsup / / dmdt = limsup / (i? n , it„)y di = / (i9,w)vdt= / / $udmdi. (142) 

nyoo JO JQ n/M Jo JO JO JfT 

Combining (|139| l- (|140|) with (|142[) and using as above [21 Prop. 1.1, p. 42], we obtain — l/fl = u a.e. in 
(0,T) x Jl and — 1/77 = ur a.e. on (0,T) x T. This concludes the proof of weak sequential stability. In 
particular, this argument could be easily adapted to provide a rigorous existence proof. 

Finally, the energy identity f|103[) can be proved simply by multiplying (|99| by 1 — 1/$ and (|101[) 
by Xt and performing standard integrations by parts. Indeed, it is easy to check that, also in the limit, 
1 — l/i? and Xt are sufficiently smooth to be used as test functions. 



4.3 Proof of Theorem 13. 4t uniqueness and regularity 

Uniqueness. Uniqueness can be easily obtained as follows: take two weak solutions (#1, Xi) and ($2, X2) 
and set (■&, X) := - d 2 , Xi - X 2 ). Then, write ([Ml) firstly for (i?i,Xi), then for (tf 2 , X 2 ), take the differ- 
ence, integrate it over (0, £), £ € (0, T], and test the result by U\ — u 2 , where Ui = — Correspondingly, 
take the difference of (| 1 1 [) and test it by X. Using the monotonicity of / and /r and proceeding as in 
[3UI Sec. 3] in order to estimate the nonlinearities coming from the quadratic terms involving X b and Ar, 
it is not difficult to achieve the desired result. The details, very similar to the proof given in [36l Sec. 3], 
are left to the reader. 

Parabolic regularization of solutions. Now we come to the proof of (|105|) - (|108|) . As before, we 
proceed by formal estimates. As noted above, also these estimates should be performed within some 
approximation scheme, whose details are omitted for simplicity. 

Second estimate. We test ([99]) by tu t = t'dt/'d 2 and add the result to the time derivative of (jlOip 
multiplied by tXt- Recalling ([80]) and using the monotonicity of / and of /r, we obtain 

^(^l|V U || 2 + ^||Vr U r||r + ^ll^llw)+t J ^ + t J^ t ul + t||VX t || 2 + t||V r Xr tt ||r 

-t J h b ut-t J^hvu T ,t < |(||Vt/|| 2 + HVrwrllr) + Q+<*)lN& 

+ t I A r (X r )X rt u r + i / A^'(X)X 2 u. (143) 
Jr Jn 
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Now, let us note that, by (|56")h we have for any e > 

t^A"(X)X t 2 M <t( £ ||X t || 2 v +c e ||X t || 2 || U ||^), ^A';(X r )X 2 ^r<t( £ ||Xr, t || 2 y r + Ce ||^llrlkr|| 2 y r ). (144) 
Moreover, 

-t /" h b u t -t J hrur,t = J tuh b + J tu r hrj + ( / u^6+ / ur^r) + ( / tu dth b + j turd t hr). 

(145) 

Now, thanks to the fact that H has zero mean value (cf. (|89p). we can modify all the integrands on the 
right-hand side by subtracting to u and Mr the quantity ran(u), as done in (|116p . For instance, the first 
pair of integrals gives 

- ^ (^J tuh b + J tu r hrj = ( J t(u - m n (u))h b + ^ t(u r - mfl(u))/i r j . (146) 

The other two pairs of integrals on the right-hand side of (|145[) arc managed similarly and then estimated 
directly in this way: 

(/ Ukb + L Urkr ) + ( / tudtkb + L tUr£>tkr ) - c !l M llvl|H||v + c<|| M ||v||a t H|| v ,. (147) 

As a consequence, taking e small enough in (I144[) , integrating (|143|) between and a generic t £ (0,T"), 
recalling (|104p and (|120[) . and using Gronwall's lemma, it is not difficult to arrive at 

^l|V^)|| 2 + ^||Vr U rWllr + ^II^WII«- / *(«(*) - m n {u{t)))h b {t) 

- j t(u r (t)-m n (u(t)))h r (t)+ f s f i? t V+ f s f rfiul 
Jr Jo Jn Jo Jr 

+ ll s||VXt||2 + ^ J ( s ll V ^r, t |lr < 0(Eo), (148) 

where, as before, the expression of Q may depend on the source H (and in particular on the additional 
regularity assumptions (| 104|) ^ . Note now that 

t(u(t) - m n (u(t)))h b (t) - J i(« r (t) - m a (u(t)))h r (t) > - ^\\Vu(t)\\ 2 - l - ||V r ur(*)llr - c > ( 149 ) 

where c = c(H) is independent of t. Hence, using again the generalized Poincare inequality of Lemma l3.11 
we finally obtain 

IMlWxv) + \\M\l~{T,T;n) + HXtllW,T;V) < Q(Vo,t-\T), (150) 

which gives the first of (|105[) and of (|107|) . Moreover, testing again (|101|) by /(X) we can write (compare 
with (Tmj)) 

ll/Wf + f T /(Xr)(/r(X r ) - 5X r + ^g- r) ' 



< jf /(x)(-X t + iX-^U) -j f(Xr)X T ,t 

< hfWW 2 + xll/(*r)|| r + <\\X\f + ||X t || 2 + ||X r , t || r ) + c(l + ||X||| 4(0) ) |M| 2 4(0) . (151) 
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Consequently, recalling (|122[) , taking the essential supremum with respect to time, and using (|5rJ[) , (|150l) , 
(|120p . and the compatibility condition (|83|) . it is not difficult to infer 

|| - AXllioo^r.H) + ||^MIU-(r,r i fl) < Q(E ,t-\T), for any r g (0,1). (152) 

Then, using again Theorem 2.7.7 and Theorem 3.1.5], we get d u X <G L°°(t, T;i?r) for any t G (0,1). 
Hence, following the lines of the argument in Section B~T1 we arrive at (compare with (|125[) - p26[0 

||-A r Xr|U-(r,T ; Hr) + ll s o,r(^r)IU-(r,T i if r ) <Q(®o,t-\T), for any re (0,1). (153) 

Relations (|152[) . f|l 53[) give ()106[) . (|108|) . respectively. Hence, to conclude the proof of Theorem 13.41 it 
remains to show the second of (|105[) and of (| 1 07[) . This essentially relics on the following result, which is 
a variant of [101 Lemma 3.3]: 



Lemma 4.1. Let (■&, 77) be a smooth solution to the problem 

1 

~~0 



d t $-Au = g b , u = --, inn, (154) 



d t i] — Ar«r = ffr ~ 9 u u, n = on V, (155) 

u T 

over the generic time interval (S,T), where we additionally assume that 

\\ u \\L3(s,T;L3/ 2 (n)) < M i ll.9b|lL 2 (s,T;L3+'(n)) + \\9r\\L*(s,T;L*+*(r)) < G, (156) 
for some (given) constants M > 0, G > and some e > 0. Moreover, let us assume that 

u(S)eL 1 (il), u r (S)eL\r), \\u(S)\\ LHn) + \\u r (S)\\ LHr) < U, (157) 
for some U > 0. Then, we have 

IMIl~((s+t,t)xQ) + IKIU»((S+T,T)xr) fCQ^Af^T- 1 ) Vre(0,T-5). (158) 

The proof is mainly based on a Moser iteration scheme with regularization and closely follows the lines 
of [401 Lemma 3.3], where the case of Neumann boundary conditions is treated. The easy adaptation is 
left to the reader. 

At this point we are in a position to conclude the proof of Theorem 13.41 Actually, assumption (|104p 
and estimate (|150p allow us to apply Lemma [4.11 with gt = hb — (Xb(X))t and gr = hr — (Ar(Xr))t on 
the generic interval (S,T) with S > 0. This provides the L°°-regularization for u and for its trace ur, as 
desired. I 

Proof of Corollary 13.51 The proof is based on the comparison principle for ODEs, similarly to [21] 
Sec. 3.2]. Given r > 0, thanks to the bounds (|105j) - (|108[) and to assumption (|5fJ|) . there exists M > 
depending on Eo and t, but otherwise independent of time, such that 

\\6X + A' b (X)u|| L o 0((T)T)xn ) + ||a*r + A^(XrHIU~((r,T)xQ) < M. (159) 
Let us now consider the following (forward in time) initial value problem 

U + Ks f(y) = M + C s , 

\y(r) = i, 
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the constants k s and C s being as in (|109[) . Thanks to (|82|l and to the standard theory of ODE's, there 
exists lu = u>(t) e (0, 1) such that y(t) < 1 — lj for all t > It. Let us then view y as a function also of the 
space variable x and subtract the first row of (|160p from (|101j) . both in the bulk and on the boundary. 
Then, test the result by (X — y) + , (-) + denoting the positive part. We have 

JL I (* - y) + 1 2 dm + / (/(*) - (* - y) + dx 

(/r(Xr) - K.f(y) + C s ) (X r - y) + dS 

(SX + \' b (X)u-M-C s )(X-y)+dx + J (6X r + \' r (X r )ur~M)(X r -y) + dS <0. (161) 

Recalling (|109[) and using in particular that k s £ (0,1], exploiting the monotonicity of /, and noting 
that (X — y) + \t= T = dm-almost everywhere in fl thanks to the fact that X takes almost everywhere 
its values in dom/ = (—1, 1), the comparison principle yields that X(t,x) < y{t) < 1 — u for all t > 2r 
and for d?n-almost all x € Q. Namely, we have the upper bound in (|110p . The lower bound is proved 
similarly. I 

Proof of Prop. 13.61 As before, we get the uniform boundedness of $ by means of a Moser iteration 
scheme. Namely, we will rely on the following lemma, whose proof can be obtained by suitably modifying 
the argument given in |401 Lemma 3.5] and is thus left to the reader. 

Lemma 4.2. Let rf) be a smooth solution of problem ()154[) - (|155j) over the generic time interval (S, T), 
where we additionally assume that 

ll^llz,-»(s,T;i,i(n,dm)) ^ M > IMIl2(s\t ; l3+=(q)) + \\gr\\ L ^s,T-L^(r)) < G, (162) 
for some (given) constants M > 0, G > and some e > 0. Moreover, let us assume that 

$(S)€L 3+ %n), v (S)eL a+ ^T), \\#(S)\\ L 3 +c( n) + \\ri(S)\\L^(n)<e, (163) 

for some (given) constant > 0. Then, for any t G (0, T — S), we have 

||l?||i-((s+r 1 r)xn) < Q{M,G,@,t- 1 ). (164) 

However, to apply the above lemma over the generic interval (S, T), S > 0, some care is needed. Actually, 
taking as above = hb — (Xb(X)) t and = hr — (Ar(X))t, (|162p is satisfied thanks to hypothesis (|104j) 
and estimate (|150p . However, it is not a priori obvious that assumption (|llip implies f|l 63[) for S > 0. 
To prove this fact, we have to provide a uniform control of the L 3+e -norm of $ over the interval (0,5). 

To this aim, we just consider the situation (more difficult, here) when — so,6 and — so.r are "regular" 
potentials and (|112p holds. In this case, we first have to prove that a uniform control of the L p -norm of 
X holds on small time intervals (this property is, of course, obvious when we have singular potentials). 
To do this, it suffices to test (|101[) by |X| P ~ 2 X, p to be chosen below, and integrate both in and on T. 
Using the monotonicity of / and fr together with assumption (|86[) . it is then easy to get, for some k > 0, 

i\W\Un, dm) + K fjn\X\^X)\ 2 + K / r |V r (|X r |^X r) | 2 

<c / \X\ p dm + c [_ (1 + |X| P )M dm (165) 
Jn Jn 
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and the latter term on the right-hand side can be controlled as follows (we just deal with the bulk 
component, the boundary one behaving similarly): 

(fi) IMI t 3 p(n) - c IMU 6 (fi) IMIlp(n) Mll3p(n) 
< ^M 2 v\\X\\ p LP(n) + a\\X\\ p L3p(n) <c a (l + \\u\\ 2 v ) \\X\\ p LP(n) + a\\ V(|X|^X)|| 2 . (166) 

Then, we first take a small enough so that the latter term on the right-hand side is controlled by the 
second term on the left-hand side of (|165[) . Subsequently, we use Gronwall's Lemma in (|165[) to obtain, 
for any p S (1, oo), 

\\X\\L~ i0 , S ;Le(TlAm)) < Q (^0. 11*0 \\ L ~ { n, dm) , S, p) . (167) 

Next, we proceed along the lines of [301 Proof of Theorem 2.7]. Namely, we test (j^)) by ?9 2+£ , both on 
the bulk and on the boundary. We then get 

^ltl { n An) - c f n \~ + h ^ +t + c J r I - A r(Xr)Xr,* + W +e 

< c(|| - X' b (X)X t + h b \\ L3+ , m + || - X'AXr)X r ,t + MU^(r)) M^l^y (168) 

whence, clearly, 

^IKIlL3+e ( n,dm) < c(||A , 6 (X)X t || i » + . ( n ) + INU»+-<n) + ||A{.(Xr)Xr,t||z»+.(r) + H^rlU»+-<r>) (169) 

and we have to estimate the right-hand side over the "small" time interval (0,5). Actually, the terms 
depending on H arc controlled thanks to assumption (|104|) . We just give an estimate for the bulk term 
depending on Xt, which is the most difficult one due to worse embeddings holding in 3D. To do this, we 
use estimate (| 167[) with p = 9 + 3e (this value is selected just for simplicity of computation) . Then, using 
assumption (|86[) and interpolation, 

\\K(X)X t \\ L ^ m <c(l + \\X\\ L * +3c{n) )\\X t \\ L ^ in) <c\\X t \\ L s^. (n) . (170) 

Thus, using that 

H^(n) C L— (Q) (171) 
and standard interpolation properties of Sobolcv spaces, we obtain 

\\X t \\^ <c||X t ||p||Xt||£* <c(\\X t \\^t^)(\\X t \\^t^)t-i, (172) 
whence, using estimate (|148p . 

l|A6(X)Xt|| £ s+.(n) <c||X t ||9±3. < ct\\X t \\ 2 v + ct- 6 ^ . (173) 

Hence, choosing e > small enough so that the latter exponent of t is strictly larger than — 1, we can 
integrate (|169p over (0, S) to obtain 

ll 1? lli«»(o,sji3+«(n I dm)) ^ <3( E 0' ll' ? o|li3+«(n idTO )) J ||Xo|| i o.(n idTO ),5'). (174) 
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Hence, (|163[) holds and we can apply Lemma |4~21 ffor S > arbitrarily small) to get (|113p . as desired. I 
Proof of Corollary 13.71 First, note that, by (|113|) . for all r > we have 

\u(x,t)\ + |« r (a;,t)| > c(r) for a.a. (x,t) G Q x (t,T), 

with c(-) : R + — > M + possibly going to when r \ 0. Hence, using (|148[) . we deduce that 

/ R(s)f ds < c{t)- 2 t~ 1 f s [ $ 2 t u 2 dxds<Q(E ,T-\T) (175) 

as well as 

y ||?7 t (s)||f-cis < c(T)-V- i y s y ?/ MdSd S <Q(E ,T-\T). (176) 

These relations correspond to (|114j) . To conclude, we need to prove that (|99[) can be interpreted in the 
strong form (|66p . Actually, this follows by reasoning as done in the "Estimate of the nonlinear terms". 
Namely, a comparison of terms in the bulk component of the heat equation gives that Au € L 2 (t, T; H). 
In turn, this provides additional regularity of d u u and, as a consequence, that ApMr S L 2 (t, T; Hr)- The 
details are left to the reader. I 



5 Long-time behavior 

In this section we prove existence and regularity properties of cj-limit sets of weak solutions to Prob- 
lem (P). Moreover, we also prove existence of the global attractor. For simplicity, we will assume that 
no external heat source is present, i.e., we will take H = 0. An asymptotically vanishing source could be 
treated as well, paying the price of technical complications. Moreover, we will restrict ourselves to the 
case of singular potentials (the case of regular potentials being in fact simpler) and, in particular, we will 
assume the reinforced compatibility condition (|109[) . Finally, we will assume here the higher summability 
property (jllip on i?o (cf. Remark 15.31 below for a motivation for this choice). 

To start with, we write the stationary problem associated to Problem (P). In order to properly state 
it, we have to notice that, testing (|99|) by 1, one gets that the quantity 

H = lt{0,X):= f ($ + \ b (X))+ f (»y + Ar(Xr)), (177) 
Jq Jr 

representing the "total mass" of the (bulk+boundary) internal energy, is conserved in the time-evolution 
of the system. Hence, once a solution trajectory evolves from some initial datum ($Oj-^o) having finite 
energy Eo, any point in the w-limit set must respect the constraint (|177p . with \x depending on the given 
initial datum. 

Moreover, it is apparent that the stationary version of (|99[) simply prescribes that u is a constant. 
Hence, we may write the steady state problem associated to Problem (P) as the following system: 

tfooe(0,+co), It*, = -l/tfoo, (178) 
-AX 00 + /(X 00 )- i 5X co ^A;(X 00 )m 00 , (179) 
— ArXoo.r + /r(Xoo,r) — SXoo.r = A r (X 00j r)itoo — <9„Xoo, (180) 

/ (tfoo + AbpCoo)) + / (^oo + Ar(X 00 , r )) =/i. (181) 
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System (|178p - (|181|) will be named as Problem (Poo,/*) m what follows. We can now state and prove our 
first result regarding existence of nonempty w-limit sets: 

Theorem 5.1. Let the assumptions of Theorem 13.41 hold and let H = 0. Moreover, let (| 109)) and 
(jllip hold. Let (i?,X) be the corresponding weak solution to Problem (P). Then, ast f~ oo, (u(t),X(t)) 
is precompact in T-L x V. Moreover, any limit point (i9oo,^oo) of (any subsequence of) ($(t),X(t)) is a 
solution of Problem (P mi)1 ), where the quantity fi in (1181[) is equal to the "mass" /i($o,Xo) of the initial 
internal energy. 

Proof of Theorem 15.11 As a first step, we need to go back to the estimates performed in Section @] 
In particular, repeating the energy estimate (|1 15|) and noting that now H = 0, we immediately obtain 

||i?-logt?|| L=0(0 +00;il(ndm)) + ||Vw|| L 2 (0:+oo;H) + ||Vru r ||L2 (0i+oo;Hr) < C, (182) 

||^IU=(0,+oo;W) + ||X||l~>(0, + oo;V) + II X IIl~((0, + oo)xQ) ^ C - ( 183 ) 

Here and below, C > denotes suitable constants possibly depending on Eo but assumed in any case to 
be independent of the time variable. Instead, we will denote by c M > the constants, also independent 
on time, that are allowed to depend on the initial data only through the conserved value /x. 
Thanks to ([TT8 | - (fTT9)) . we also get 

ll«IU»(t,tfl;V) < C. (184) 

In particular, (fl82 ]) -([l83 |) imply that 

Vt>0, 3Se [t,t + l] : ||u(S)||v + ||X t (S)ll*<C- (185) 

Then, we can go back to the "Second estimate" of Subsec. 14.31 Actually, we test (|99| by u t and the time 
derivative of (|101|) by X* (hence, we do not need here the weight t). This leads to the analogue of (|143p . 
namely, 

^(^l|V U f + ^||V r ur|lr + l\\X t f H ) + JjW + J^Wv + l|VX t || 2 + ||VrXr,t||£ 

< / K{X)X 2 t u+ f A?(Xr)(X r . t ) 2 «r. (186) 

Then, estimating the above right-hand side as in (|144[) . integrating over (5, 5 + 2), S as in (|185p . and 
using the generalized Poincare inequality (|85p . we readily get 

||X t ||i=o( 1)+00 . W ) + ||u||l~(i,+oo ; v) < C, (187) 
\\(\ogd) t \\ L 2^ t+1 . u) + ||X t || i2(tit+1;V) <C, Vt e [1, +oo). (188) 

Going back to (|186p . and having now (|187p at our disposal, we note that the V-norm of u on the right- 
hand side of (|144p is controlled uniformly in time. Hence, we can integrate (|186p over the whole (1, +oo), 
so that (|187p is improved to 

ll^ll^(l,+oo;V) + ||Oog0)t|Ua(l >+ oo;W) < C. (189) 

Next, we can repeat estimate (|151[) over (5, 5 + 2), obtaining 

||/(X)IU«( ll+00i J?) + ll/r(Xr)|Uoo (1 , +oo;Hr) < C. (190) 
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This allows us, as before, to get also 



IMU°°(l,+oo;if 2 (£2)) + ll^r||L°°(l. + oo;// 2 (r)) < C. (191) 

Next, (|182[) - (|183[) and a comparison of terms in (IM1) allow us to see that 

||0t||i»(O,+oo;V) < C. (192) 

Thanks also to the Aubin-Lions lemma, the first of (|187[) . (|191|) . and the compact embedding W C V, 
immediately give the precompactness of X(t) in V. 

Applying Lemma |4~T1 over the generic interval (S, S + 2), we also infer that 



ll«lli~ ((li+oo)xn) < C, (193) 

which, combined with (|189[) , gives 

IMIl 2 (1, + oc;«) <C. (194) 

Thanks to (|193j) . we can repeat the comparison argument of Corollary 13.51 In particular, we get the 
separation property uniformly in time: 

\X(t, x)\ < 1 — u> dm-almost everywhere in f2 for a.e. t > 1, (195) 

with uj 6 (0, 1) independent of time. 

By (|194[) . (|187|) and the Aubin-Lions lemma we deduce that u(t) is precompact in H. With the 
precompactness of (u(t),X(t)) and the above estimates at our disposal, we can now prove that the 
trajectory (#(£), u(f)) admits a nonempty w-limit set. 

Namely, we take {t n } to be a diverging sequences of time and we consider Problem (P) over the time 
interval (0, 1) and with initial datum ($(£„), X(t n )). Let us call Problem (P„) this problem. Hence, once 
is a weak solution to (P), (#„(*), X n (t)) := (i?(i n + t), X(t n + *)), i 6(0, 1), solves (P„). 

Thanks to the previous estimates and to the Aubin-Lions lemma, there exist limit functions u, X 
defined over (0, 1) x such that, say, 

u n strongly in C°([0, 1]; F 1_e (0) x H 1 ' 5 ^)), (196) 
X n -> X strongly in C°([0, 1]; V), (197) 

for all e > 0. As usual, all convergences arc to be intended up to the extraction of (nonrclabcllcd) 
subsequences of n /* oo. Moreover, thanks to (|1 89[) and (|194[) . both u and X are constant in time; hence, 
they coincide with the limit of (the extracted subsequence of) u(t n ) and X(t n ), respectively. Moreover, 
thanks to (|182j) . u is also constant in space and, thanks to (|193[) . we have |u| < C for some C > possibly 
depending on E . 

Hence, passing to the limit in (|101|) . we get that (u, X) solve (|1T0|) - (|180[) . which can be written in the 
strong form thanks to the same considerations on regularity as those made for the evolutionary system. 
In particular, the limits of the terms f(X n ) and f(X n .r) can be identified as before by monotonicity 
methods. 

To conclude the proof, we have to take the limit in the heat equation (|99| in order to recover (|178j) 
and (|181|) . To do this, we recall that, for all t <E [0, 1], we have 

Mt) = Vn(t) = — ( * . (198) 

Unit) {Un(t))r 
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almost everywhere in fi and, respectively, on T. Then, by (|196[) , we have 

n (t) -)• -i, ffc,(t) -=31, (199) 

u u 

almost everywhere in f2 and, respectively, on T. Let us now prove that, actually, u cannot be 0. Indeed, 
by (|177p , we have, for almost all i £ [0, 1], 

/ ■& n (t)dm<\fi\ + \n\ max |A&(r)| + |r| max |A r (r)| < c„. (200) 
J?T re[-i,i] re[-i,i] 

Consequently, applying Jensen's inequality, we deduce that there exists > depending only on /i, 
and Ar such that |u| > c M . To characterize the limit of we need the following lemma: 



Lemma 5.2. Under the assumptions of Theorem 15. 1[ if we have in addition 

(i?o, %) G L p (n) x L p (r), for some p > 3, (201) 

then it follows that 

\\(m,v(t))hnn)xLnr) < Op(Eo, ||<?o||i»(ii), ho|Up(r)), Vt > 0, (202) 
with Q p independent of t S [0, +oo). 

Proof. For simplicity, we just give the proof in the case when p = 3, the general case following by 
repeating the procedure and applying a simple bootstrap argument. Firstly, let us observe that, by (|174j) 
with e = 0, we have 

HW*),»/(t))IU»(n)x£»(r) <Q(Eo,||^o||L»(n),||»»IU»(r)) Vie [0,1]. (203) 

Hence, we can test (|99p by -d 2 . Using that now X is uniformly bounded, we obtain (compare with (| 1691) 1. 
for some k > (which may change from line to line in the computations below), 

^11^1^3^^^ + ^HV^ 1 / 2 !! 2 + ^ll^r^ 1 / 2 !!^ < c ^_|^ £ ^ 2 | d^,. (204) 

Actually we can notice that, if (|201[) holds, then the solution is regular enough so that r\ is in fact the 
trace of -d in this case. 

Noting that, due to (fT77|) . 

Htf 1/2 (*)llW,cM = j_m&™ < c M Vt G [0,+oc), (205) 
adding (|205[) to (|204p . and using standard Sobolev's embeddings, we infer 

Tull^llL^TLdm) + K H^!lL 3 (adrn) - C ll X *lll,3(n,dm)ll^ 2 |lz, 2 (f2,dm)ll^ 1/ ' 2 |l£ fi (fT,dm) + C M 

< 4M\l m>Am) \\nl^ Am) + f H*llL»(n,d») + V (206) 

Thus, setting 

y(t) := max{l, ||0(t) || 2 L3( ^ dm) }, m(t) := ||X t (t)||^ ( ^ dm) , (207) 
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it is clear that (|206[) can be interpreted as the differential inequality 

y'(t) + K < cm(t)y(t) + c^t)" 1 ' 2 , (208) 

so that the comparison principle, together with (|203[) . (|189|) and Gronwall's lemma, readily imply that 
y is uniformly bounded for t G [l,+oo), as desired. Indeed, in the set of times such that y(t) < Ac^/k 2 
there is nothing to prove, while if y > Ac^/k 2 (|208[) takes the form y' < amy and we can use Gronwall's 
lemma since m is globally summable thanks to (|189[) . I 

Remark 5.3. It is worth noticing that the property stated in the Lemma is likely to be false if p G (1, 3). 
Indeed, in that case the gradient terms in (|204j) seem to give no help and, consequently, the differential 
inequality corresponding to (|208|) loses its dissipative character (in other words, it provides a control on 
the _L p -norm of $ which is not uniform in time) . This is the reason which led us to assume (jllip in this 
section. On the other hand, thanks to the conservation of \i (cf. (|177[1 ), we have in any case a uniform in 
time L 1 -control on However, this seems not enough in order to characterize properly the limit of $ n . 
Indeed, it may happen that 

l?n = #n,l + $n,2, so that u n = — — , 

where is a "good" function which converges, say, uniformly, to the constant —1/u. Instead, the 
sequence {$71,2} is also bounded in L 1 but it may converge to some singular measure (e.g., a Dirac mass, 
meaning that the support of $ ni 2 becomes small and i?„ ; 2 becomes large in that set). However, when one 
computes u n , the contribution of $ nj 2 is negligible for large n because $ nj 2 influences the value of u n only 
in a set that has asymptotically measure 0. Hence, the limit of d n is in this case a measure whose regular 
(absolutely continuous) part coincides with —1/u, but which may also have a singular component. 

Thanks to Lemma 15721 using pointwise convergence (cf. (|199[) ) and Lebesgue's theorem, we obtain 

K^-^t strongly in L p (0, 1; L P (H; dm)) VpG[l,3). (209) 
u 

In particular, we get the second (|178[) (|18ip . which concludes the proof of Theorem 15. 11 I 

Finally, we come to the problem of the existence of the global attractor for Problem (P). To address this 
issue, we first introduce the natural phase space for the dynamical system associated to Problem (P), by 
setting 

{(& > dm-a.e. in fi, logti G L^U; dm) 1 
(0, X) G (L 3+e (r>) x L 3+e (r)) xV: F(X) G L 1 (f2), F r (X r ) G L 1 ^), \ , (210) 

[f Q {# + \ b (X))+f r (r l + \r(Xr))=li J 

where e > is given and the conserved value n is also prescribed. It is not difficult to prove (see, e.g., 
[3SJ for details) that the above space, endowed with the distance 

dirt*((0i,Xi),(0 a ,X 3 )) := \Wi-M L ^(n,dm) + ll lo g^- lo g^|| il(TT , dm ) 

+ ||Xi - X a || v + ||F(Xi) - F(X 3 )|| £ i (n) + l^r^i.r) - F r (X 2X )\\ L1(r) (211) 

acquires a complete metric structure. 

We will name Problem (P M ) the version of Problem (P) where the "mass" \i of the initial internal 
energy is assigned. 
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Remark 5.4. In order for the phase space X^ not to be empty, we implicitly assume that 

H > 101 min A 6 (r) + |r| min A r (r). (212) 

re [-1,1] rG[-l,l] 

Actually, if (|212[) is not satisfied, then from (|177p we get that # needs to be nonpositive, and the problem 
becomes inconsistent. However, we will see just below that a stronger condition may be needed. 

Theorem 5.5. Let the assumptions of Theorem 1 5 . 1 1 hold. Let also assume that, either 

H > 101 max A 6 (r) + |r| max A r (r), (213) 
re[-l,l] re[-l,l] 

or 

liminf A'(r) signr > 0. (214) 

Then, the dynamical process associated to Problem (P M ) admits the global attractor A^. Moreover, there 
exist constants K > 0, u) G (0, 1) and a G (0, 1) depending only on u. such that, for any (■#, X) G .4^, 

Mv + IMIw (215) 
a < i9 < a -1 dm-ahnost everywhere in O, (216) 
— l+w<X<l — uj dm-almost everywhere in O. (217) 

Proof. Due to the singular and degenerate character of equation (|99p . a direct proof of a dissipative 
estimate (e.g., of the existence of an absorbing set) appears to be out of reach. On the other hand, as 
in |39| , we can take advantage of the fact that our system admits a coercive Liapounov functional (the 
energy 8). Thanks to this property, the existence of the global attractor follows by proving the following 
conditions: 

(a) Solution trajectories are precompact with respect to the metrics dist^ of the phase space X^\ 

(b) The set of stationary states is bounded in the energy space independently of the magnitude of the 
initial data. 

Hence, let us start with the proof of (a). By Theorem 15.11 it is clear that for any weak solution there 
exists C depending on Eo and on the L 3+e -norm of $0 such that 

||u(t)|| v +||X(t)|| w <C Vte[l,+oo). (218) 

Moreover, thanks also to Lemma r5.2[ we have 

< a < < a" 1 dm -a.e. in O and for all t G [1, +00), (219) 

with a G (0, 1) depending on the same quantities as the above constant C . Coupling (|218[) with (12191) . 
we readily get 

||0(t)|| v < C* Vt G [1,+oc). (220) 

In addition to that, the uniform separation property (|195[) holds. Hence, both ■& and X are eventually 
separated from the singularities. Using (|218[) - (|220[) and (|195[) it is a standard matter to prove that the 
trajectory is precompact with respect to the metric structure of X^. 
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To conclude the proof, let us demonstrate (b). This is more delicate, since we need to prove that 
the set of stationary states (i.e., of the solutions ($ooj ^oo) to Problem (Poo, /J) is bounded in a way that 
depends only on /i. In other words, all elements of the w-limit of a given solution trajectory are bounded 
in a way that depends on the initial datum only through the conserved quantity /i. Actually, we have 
already seen (see (|200p ) that \"doo \ < c M . To prove that the same property holds for Uqo, we first consider 
the situation when f|2 13[) holds. In this case, combining (|213[) with (|18ip . \uco\ < c n follows immediately. 

Hence, the terms X'^X^Uoo and Ap(Xoo,r)woo,r on the right-hand sides of (|179p and, respectively, 
(|180p arc uniformly bounded in a way depending only on fi. Then, by the maximum principle it immedi- 
ately follows that there exists w M e (0, 1) such that — 1 + uj^ < Xoo(t, x) < 1 — w M dm-almost everywhere 
in O. It is immediate to check that the same property holds if (|2 14[) is satisfied in place of (|213|) . 

Standard elliptic regularity estimates applied to (|179[) - (|180p then give, say, 

||Xoo||w<c M (221) 

for every stationary solution ($00,^00) (actually, much more is true since X x is separated from the 
singularities and is, consequently, a classical solution to the elliptic system (|179p - (|180p ). Hence, (b) holds 
true and the existence of the attractor with the properties (|215|l - (|217|) follows from classical results (see, 
e.g., M)- ■ 
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